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Abstract 

The anomalous Hall effect in two-dimensional ferromagnets is discussed to 
be the physical realization of the parity anomaly in (2+l)D, and the band 
crossing points behave as the topological singularity in the Brillouin zone. 
This appears as the sharp peaks and the sign changes of the transverse con- 
ductance function of the Fermi energy and/or the magnetization. 
The relevance to the experiments including the three dimensional systems is 
also discussed. 
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Anomalous Hall effect (AHE) in ferromagnets is a phenomenon that there 
appears the contribution to the transverse resistivity pn proportional to the 
magnetization M in addition to the usual Hall effect |[J: 

p H = R B + At:R s M (1) 

where B is the magnetic field, Rq is the ordinary Hall constant, and R s is 
the anomalous Hall coefficient. It has been recognized that this AHE comes 
from the spin-orbit interaction A because some connection between the spin 
and orbital motion is needed [| |L 0, BJ. Karplus and Luttinger (KL) 
[0] was the first to study theoretically the AHE based on the band theory, 
and recognized that this effect is due to the inter-band matrix element of 
the current operators. An important point is that the eigenstates of the 
Hamiltonian are still the Bloch waves \nk) (n: band index, k: crystal mo- 
mentum in the 1st Brillouin zone) even when the spin-orbit interaction is 
present, namely the translational symmetry remains. Therefore a xy (u) for 
noninteracting electrons is given by [|7| 



a xy {u) = i 



EE 



f(e m (k))-f{e n {k)) 



^ (nk\ J y \mk)(mk\ J x \nk) ^ 
lo + ir] + e m (k) - e n (k) 

where J M is the current operator, f(x) is the Fermi distribution function, and 
we take the natural unit, i.e. h — c — 1. This is the general expression, and 
KL result can be obtained by expanding the above Kubo formula to the first 
order in A. On the other hand, Smit stressed the essential relevance of the 
impurity scattering, which disturbs the translational symmetry. Luttinger 
||] reconsidered this problem by taking into account the impurity scattering 
in a systematic way, and essentially found that o xy = (T KL + a L where <jkl is 
the original result of KL || while the second term is the new term which is 
proportional to the inverse of the impurity scattering potential, is strongly 
dependent on the explicit form of the impurity scattering, and vanishes for 
the short range potential |J . However most of the later theories regards the 
modification of the scattering events due to the spin-orbit interaction as the 
origin of AHE [|. 

Very recently, it has been experimentally recognized that there are two 
distinct types of behaviors of AHE || £J |ID[ [□], [L2|, [L3| . In the conventional 
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one, a xy decreases as the temperature T is lowered, and the AHE is attributed 
to the scattering events influenced by the spin-orbit interaction (type I). In 
the other case, a xy increases as T is lowered and remains finite in the limit 
T = (type II). In the type II systems, various anomalous behaviors of 
AHE are observed. For example, pu in thin films of SrRuOa shows the 
change of magnitude of 10~ 9 f2m with the sign inversion [lOj as a function of 



the temperature T below the ferromagnetic transition temperature T c . The 
change of a xy per layer is estimated as 0(e 2 /h) by using the data of the 
longitudinal resistivity p ~ 10~ 7 f2m and the averaged lattice parameter a ~ 
4 x 10~ 10 m. (We appropriately recover Planck constant h.) Another example 
is the pyrochlore ferromagnet (Snii-zCa^Mx^Oy showing the strong doping 
dependence of AHE, where hole doping of the order of x = 0.1 causes the 



factor of ~8 change of pn in contrast to the almost unchanged p xx [[L2 
Similar phenomenon is observed also for the temperature dependence of pn 
in below T = IK (< T c ^ 100K) in Nd 2 Mo 2 7 @. This huge sensitivity on 
Ep and/or T can not be understood in terms of the conventional theories on 
the AHE nor the change of Rq below T c . The purpose of this paper is to give 
a theoretical explanation for the distinction of type I and type II systems 
and the anomalous behaviors of a xy in type II systems. 

On the other hand, the topological nature of a xy in 2D has been exten- 
sively studied in the context of quantum Hall effect under external magnetic 
field ||14|| . The Kubo formula eq.(0) for uo = can be rewritten as ]14 



n.k 



(nk\J y \mk) (mk\J x \nk) — (J x <-> J y ) 
(k) - e m {k)] 2 

= e 2 ^/(£#))[v t xl n (i| (3) 

n,k 

where A n (k) = —i (nk\Vk\nk). When the Fermi energy is in a gap and 
T = OK, this integral is transformed to the "vorticity" corresponding to the 
vector potential A n (k), and is an integer called Chern number [times (27r) -1 ]. 

From the viewpoint of the field theory, the appearance of a xy in 2D 
is closely related to the parity anomaly in (2+l)D Dirac fermions, where 
the Chern-Simons term is generated with the sign depending on that of the 
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fermion mass In the Bloch wave case, this Dirac fermion corresponds to 
the band crossing [16], and several models have been proposed for the con- 
densed matter realization of the parity anomaly [T7|, |I8| . Haldane [TPS] recog- 



nized that the quantum Hall effect can be realized as the result of the mag- 
netic ordering and without external magnetic field in 2D. It should be noted 
here that parity transformation in (2+l)D is r = (x, y) — > f = (—x,y), and 
this symmetry operation changes the sign of the Chern-Simons term. There- 
fore the effective Hamiltonian should break both P(2+ip- and T-symmetries 
(at least infinitesimally) . 

Although the physical realization of Haldane's model i.e., the com- 
plex transfer integrals between next-nearest-neighbor sites on the honeycomb 
lattice, is not so easy, more realistic model has been proposed in connection to 
the the spin chirality mechanism of AHE [pl| . The spin chirality Si ■ Sj x Sk is 
basically the solid angle subtended by the non-coplanar spin configurations, 
and its relevance to the AHE has been investigated in manganites at finite 



temperature || and in pyrochlore compounds || [IT] at zero or low temper- 
atures. This spin chirality acts as the effective magnetic field [p^] , and it 
has been found that the geometry of the lattice is crucial for the finite a xy 



in terms of this mechanism [19]. Roughly speaking, the positive and nega- 
tive fluxes in the unit cell should be inequivalent to avoid the cancellation. 
The multiband effect due to the existence of plural number of atoms in a 
unit cell is essential, and each band is characterized by the Chern number 



14j| . Therefore the two conditions, i.e., non-coplanar spins and non-trivial 
geometrical structure, makes this mechanism rather special. 

In this letter we explore the connection between the AHE and the par- 
ity anomaly in 2D ferromagnets, and found the followings. (i) The parity 
anomaly and the resultant nonzero a xy is not the special case but is the 
generic phenomenon occurring in P(2+i)d~ and T-broken magnetic ordered 
state. Our model is the conventional one with the spin-orbit interaction A 
and the collinear spins on the simple square lattice (perovskite structure), 
(ii) The band crossings occur very often as we change the spin-orbit coupling 
A and/or the magnetization m z . Even though the Fermi energy is in general 
not in a gap, these band crossings (Dirac fermions) behave as the topological 
singularities and appear as the sharp peaks and/or the sign change of a xy as 
a function of the Fermi energy or the magnetization m z . 

We consider a model with triply degenerate t 2g orbitals with spin-orbit 
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interaction on the 2D square lattice. 



H = - ^t Q [cl ya (i)c xya (i + x) + cl ya (i)c xyu (i + y) 

+c\ za ii) c yzaii + y) + e\ xa {i)c zxa {i + x) + h.c] 
+ Yl ±t d4 y *( i ) c zx*(i±y) + c xy(T (i)c yza (i±x)+h.c.} 

+H so + H el _ el (4) 

where the transfer integrals are determined by the oxygen p-orbitals be- 
tween the two transition metal ions, to is nonzero even for the perfect per- 
ovskite structure, while ti becomes nonzero due to the shift of the oxygen 
atoms out of plane in the positive (negative) z-direction. This will cause 
the Dzyaloshinsky-Moriya interaction |2l| when the spin-orbit interaction 
Hso and the transfer integrals are combined, because the space inversion 
symmetry at the middle of the bond is broken. Our definition of P(2+i)d 
here is for the effective 2D Hamiltonian, and does not include the trans- 
formation on the orbital and spin degrees of freedom. The ti-term breaks 
this P(2+i)D-symmetry. If the shift of the oxygen atom is alternating, the 
sign ± in the t\ term is missing and P(2+i)D-symmetry is unbroken. In 
this case, the topological effect described below is missing, because the 
gauge flux density b n (k) = [V* X A n (k)] z in eq.(|3]) has the symmetry as 
b n (—k x ,k y ) = —b n (k x , k y ). (Note that b n {k) itself is not zero.) We approxi- 
mate Hso by taking the matrix elements of the angular momentum I in the 
space of t 2g orbitals as 

H so = ^EE c L(04/3 • WC^W (5) 

where r are Pauli matrices and a, (3 = xy, yz, zx. Although Hso is the on-site 
interaction, it can produce the effective gauge flux combined with the transfer 
integrals. The electron-electron interaction i? e i-ei induces the ferromagnetic 
ordering and its details are not relevant to the present study. We take the 
mean field approximation as 

#el-el = ~J m z O'CLWCWW, ( 6 ) 
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where the ferromagnetic moment is along the z-direction. Now the 6x6 
Hamiltonian matrix for each fc-point can be easily diagonalized, and 6 bands 
are obtains. Figure |T] shows the energy dispersion obtained for a set of 
parameters given in the caption. The near-degeneracies of the two bands are 
found at k = [0,0] and [±vr/2, ±vr/2]. 

Fig.||(a) shows the typical example of function of Ep for the 

same set of parameters as in Fig.[l], while Fig.|](b) shows the variation as 
a function of the magnetization m z with the fixed Ep = 0. Both of them 
vary rather rapidly with the magnitude of e 2 /h and have sharp peaks. To 
analyze this behavior, we first define the Chern number Ch n for each band 
n as Ch n = i st BZ ^[^k x A n (k)] z . For the parameter set of Fig.[l], the 
ChnS are (-1, -2, 3, -4, 5, -1) from the lower to higher energy band. Namely 
the spin-orbit interaction gives the non-trivial topological structure to the 
Bloch bands in the ferromagnets as in the case of quantum Hall system. 

The next question is how the topological number appears in a xy . In the 
absence of the gap, one has to consider the distribution of [V& x A n (k)] z in 
the /c-space. In Fig.|3| shown [V*. x A n (k)] z again for the parameter set in 
FigsJI] and 0(a), where sharp peaks near k = [0, 0] and k = [±7r/2, ±7r/2] are 
clearly seen. This behavior can be understood as follows. When the band 
degeneracy occurs at k = k for u = Um z = u c or A = A c , the electrons near 
this point can be described by the generalized Dirac fermions as 



H= I -0- 2 ^(k)H(k)i;(k), «- = J-(£-£ ), 



H(k) 



V(k)+mv 2 vkp,{n x — iKy) p 
vkn(K x + i^yY V(k) — mv 2 



(7) 



where if)(k) is the two component spinor corresponding to the two crossing 
or touching bands, v and kp, are positive constants, and V{k) is a function of 
k. The usual Dirac fermion corresponds to the case p = 1. By diagonalizing 



H(k), we get the dispersion around k as e±{k) = ±J (vkp,K p ) 2 + (mv 2 ) 2 + 

V(k), where the sign +/— represents the upper/lower band. The current 
operator is given by J = J d 2 k /(2ti) 2 ^(k)VkH{k)ip{k). The mass m is a 
continuous function of (A, u), and can changes sign across the critical lines in 
(A,w)-plane. There the transfer of Ch n occurs between these two bands [p^ 



Near this criticality, i.e., when the mass \m\ is small, the gauge flux density 
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b±(k) = [Vfc x A±(k)] z is peaked near the degeneracy point as 

b / A (v { (±, fc| J y |=F, fc)(=F, fc|Js|±, fc) ~ (js Jy) 

k ± (£)-%(£)] 2 

(pvK, p ~ l ) 2 mv 2 



2[(vk D nP) 2 + (mv z 



(8) 



The flux density b±(k) contributes ±| sgn(m) to the Ch n of each band. 

The transfer of Ch n has been studied in connection with the usual Dirac 
fermion, i.e., p — 1, mainly in the context of the integer quantum Hall 
effect ||23|1. However, in the present model, we have actually found also its 



generalizations given by eq.(|7|). For example, the type of p = 2 is seen with 
the parameter set (ti = 0.5i , A ~ 0, Um z = 5t ), and the transfer of Ch n 
at each degeneracy point is ±2 when A is changed across zero. Each band 
touches the next band with the energy dispersion E(k) oc k 2 . Another remark 
here is that the band crossing and the transfer of Ch n occurs frequently. It 
occurs at least 6 times for the parameter set (t\ = OMo, A = 0.4t ) when 
Um z is changed from to 5to- It occurs more often if we additionally change 
A. 

Some remarks are now in order. One is on the effect of disorder, which 
has been neglected thus far. It is well known that the extended states in 2D 
is unstable against the infinitesimal disorder in the absence of magnetic field. 
On the other hand, the topological Ch n should be quantized and stable as 
far as the impurity scattering strength is weak enough compared with the 
gap of the Dirac fermions. This is the basis of the existence of the extended 



states and plateau of a xy in the integer QHE, and should work here also [^4 
Therefore we expect that the random potential realizes the quantum Hall 
effect in 2D ferromagnets. 

Second remark is on the three dimensional case. The crossing of bands 
also plays an important role here, but the singularity is more mild. According 
to the analysis of Dirac fermions in 3D, the peak of a xy near the critical point 
is of the order of \m\ and the width is again \m\. The magnitude of a xy is 
of the order of \m\/t Q in this case, but it varies over the energy scale of to 
except near the singularity. Near the singularity, a xy changes rapidly with 
the energy scale of |ra| while the magnitude remains of the order of |m|/i - 

Now the experimental relevance of the present results is briefly discussed. 
The most significant consequence of the topological nature of AHE is its 
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sensitivity to the Fermi energy Ep and/or the magnetization, in contrast to 



the diagonal conductivity a xx [25]. This naturally explains the anomalous 



features of AHE in thin films of SrRuC>3 |TIJ and the single crystal of py- 
rochlore Mo oxides [|12|, |13[ mentioned on the introduction. Especially it has 
been known that the topological mechanism for AHE is relevant in the latter 



material [11]. The 2D system is even more interesting, where the variation in 



a xy of the order of e /h, and with the impurity scattering the quantum Hall 
effect can be also expected. It is noted here that the conventional formula 
eq. (|^) does not work in this case, and a xy should be plotted instead. In con- 
clusions, we have studied the anomalous Hall effect from the viewpoint of the 
topology and parity/chiral anomaly. When both the parity and time-reversal 
symmetries are broken, the AHE survives even at T = where the period- 
icity is recovered ( type II). In this case, the band crossing behaves as the 
topological singularity causing the anomalous behavior of (7 X y £L function 
of the temperature T, the magnetization m z , and/or the Fermi energy Ep. 
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Figure 2: (a) Transverse conductivity function of the Fermi energy 

Ep with the parameter set t\ = 0.5to, A = 0.4t , Um z = 0.95to- The sharp 
peak near — 2t is due to the near-degeneracy point at k = [0, 0] between the 
4th and 5th bands (see Fig.|IJ), and, in other words, due to the sharp peak in 
the gauge flux density in /c-space (see Fig.||). (b) Transverse conductivity a xy 
as a function of the magnetization m z [times U] with t\ = 0.5to, A = 0.4t , 
Ep = 0. 
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Figure 3: Gauge flux density in fc-space of the 5th band with the parameter 
sets t\ = 0.5to, A = 0.4i , Um z = 0.95to for the upper, and Um z = 1.05to 
for the lower. The transfer of Ch n occurs between the 4th and 5th bands 
at (Um z ) c ~ l.Oio- 5 peaks are found in both parameter sets. However, the 
transfer occurs only at k = [0, 0] in this case. 
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